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Abstract 

We formulate the basic properties of q-vertex operators in the con- 
text of the Andrews-Baxter-Forrester (ABF) series, as an example of 
face-interaction models, derive the (/-difference equations satisfied by their 
correlation functions, and establish their connection with representation 
theory. We also discuss the g-difference equations of the Kashiwara-Miwa 
(KM) series, as an example of edge-interaction models. 

Next, the Ising model-the simplest special case of both ABF and KM 
series-is studied in more detail using the Jordan- Wigner fermions. In 
particular, all matrix elements of vertex operators are calculated. 



1 Introduction 

In this paper, we continue our study of the role of vertex operators and differ- 
ence equations in computing correlation functions in off-critical solvable lattice 
models. We start with a review of the conformal field theory (CFT) approach 
to computing critical correlation functions in WZW-type models: models with 
an affine Lie algebra symmetry, then compare it with our approach in solvable 
lattice models of the vertex type. 

At criticality 

In CFT, the correlation functions on the projective plane P 1 : 

(A |^' yl (Ci)---$xr 1,y "(C»)|A„) (1) 
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are given in terms of vertex operators acting as intertwiners in the following 
sense: 

^' Vi (0--V(\ j )^V(\ j - 1 )®Vi (2) 

where V(A) is a highest weight representation with highest weight A and highest 
weight vector |A), and is a finite dimensional representation with a spectral 
parameter £. For simplicity, we take all V J 's to be isomorphic to the same finite 
dimensional representation V. 

The correlation functions satisfy the KZ equation, a system of differential 
equations with respect to the variables £i,...,£ n . The solutions to the KZ 
equation have a branch point at Q — and their monodromy is dictated by 
the commutation relations of the vertex operators 

^• vl (Ci)< 4 ' v2 (C 2 ) = T,<' v2 (^ vl (^ w (x 2 xX (3) 

which is understood in terms of an appropriate analytic continuation, and the 
connection coefficients W are constants. 



Off criticality 

Away from criticality, one expects appropriate generalizations, or deformations 
of the structures that appear at criticality. In Jl2[ , g-deformations of the vertex 
operators and the KZ equation were obtained. There is a number of major 
differences between the critical, or classical (i.e. q = 1 ) and q-deformed cases. 
To start with, the g-deformed KZ equation is a system of difference rather than 
differential equations. Furthermore, the commutation relation changes to 

i? yl ' y2 (Ci/c 2 )<^ 1 (Ci)^r v2 (c 2 ) 

:< y \c2)*t vl (Ci)w(^ x * C1/C2V (4) 



thus, the i?-matrix appears on the left hand side, and the connection coeffi- 
cients W now depend on the spectral parameters. In the g-deformed case, the 
multivaluedness reduces to an infinite set of poles at Cj/Cfc = q n for certain 
integers n^O. 

Another crucial difference is that, in the g-deformed case, it matters whether 
the finite dimensional part Vj! appears as the right or left component in the 
action of the vertex operator: there are two different vertex operators, depending 
on the position of in the action. In (|, Q, they are distinguished as type 

I (V^ is to the right) and type II (V^ is to the left) vertex operator. In this 
language, (2) shows the action of a type I vertex operator. The operators that 



2 



appear in (4) are also of type I. They are used to compute the correlation 
functions in the form 

tT V{Xo)q - 2 ^ V \C 1 )---^ Vn (Cn), (5) 

where p is a certain grading operator, and Ao = A„. The above form of the 
correlation functions suggests a connection to conformal field theory defined 
on the annulus. Recently, an interesting paper appeared in this direction po| , 
where the classical and twisted versions of the trace function were found to 
satisfy the classical and elliptic versions of the KZ equation, respectively. On 
the other hand, the Type II vertex operators serve as creation and annihilation 
operators that diagonalize the Hamiltonian. 



Previous results 

In our previous work on the off-critical correlation functions we started with 
the XXZ spin chain, then considered the 8-vertex model. In the case of the 
XXZ spin chain, we followed the representation theoretic approach, since that 
was possible. The trace correlations that we obtained, as in (5), also satisfy a 
version of the g-deformed KZ equation (see (5) in |L5j) although the derivation 
in this case is different from that in the case of correlations that are matrix 
elements as in (1): it relies on the cyclic property of the trace and the following 
commutation relations of the vertex operators, which are valid in the principal 
grading for V z : 

= *aJ v (C7C)C, ( 6 ) 
i? vl ^ 2 (Ci/C 2 )<' yl (Ci)^' y2 (C 2 ) = ^' y2 (C 2 )^' vl (Ci). (7) 

Compare (4) with (7) : The reason that the summation on the right hand side 
disappears is that the level of the highest weights and the finite dimensional 
representation V are such that that the highest weight A in 

^(0--V(X')^V(X)^V c , (8) 

is uniquely determined from A'. Specifically, in the case of the XXZ spin 
chain, we have: level = 1, and diml^ = 2 for sl(2); refer to the definition of 
perfect crystals in |jl8| . This trivializes the connection coefficients W, that would 
otherwise appear on the right hand side. 

The above simplification was used in the derivation of the (/-difference equa- 



tions in the the 8-vertex model (see (6) in JTaj). In 15 , an argument without 
explicit use of representation theory, for deriving similar equations in the con- 
text of the 8-vertex model was proposed. The idea in this case is an extension 
of Baxter's corner transfer matrix (CTM). We replace V(A) with the space THS 1 ' 
spanned by the eigenvectors of the CTM with the specific, say i-th boundary 
condition, q is regarded as the crossing parameter, and £~ p as the CTM. The 
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basic idea is to interpret the vertex operator as a graphical insertion of a half 
line. By doing that, the commutation relations (6) and (7) follow from a simple 
graphical argument, and as a corollary the elliptic difference equation is shown 
to be valid for the correlation functions of the eight- vertex model. A mathemat- 
ical setting for the eight-vertex model, in which the vertex operators are given 
a mathematical basis, in terms of e.g. representation theory, is still unavailable. 

Outline of results 

Now, we turn to the subject and results of the present paper. We wish to pur- 
sue both the graphical and the representation theory approaches to the vertex 
operators and difference equations, but this time in the context of the ABF 
models. 

In Section 1, we go through the graphical procedure followed in the case of 
the 8-vertex model, in order to obtain the vertex operators and the (/-difference 
equations of the ABF models. The commutation relation reads as 



where the indices i and I label the different sectors in the space of states, and 
the heights on the lattice, respectively as will be explained later. Next, using 
the results of |L6), we give a representation theoretic realization of the vertex 
operators satisfying the above commutation relations. 

At this point, let us emphasize the following: The graphical derivation of the 
difference equations applies to other solvable models, even when the represen- 
tation theoretic setting is not available. Examples are edge interaction models, 
such as the Kashiwara-Miwa and Chiral Potts models. We give a brief deriva- 
tion of the difference equations of the Kashiwara-Miwa model. The treatment 
is basically a special case of that of the ABF models. 

In the case of the XXZ spin chain, we have, not only the representation 
theoretic realization of the vertex operators, but also a more powerful realization 
in terms of free bosons. In [O], an integral representation was obtained for the 
correlation functions using bosonization. In the ABF case, we do not know of a 
bosonization, or any other tool, that would provide us with an integral formula 
for the correlation functions. However, in the case of the Ising model, which 
is an intersection of all the models discussed above, with the exception of the 
XXZ spin chain, an alternative tool is available: the model can be formulated 
in terms of free fcrmions. 

The diagonalization of the Ising CTM was given in using Jordan- 

Wigner fermions. In this case, the role of type II vertex operators is basically 
played by the diagonalized fermions. In Section 2, we use this formalism to 



$ (M+1) (Ci)h,^ (4+M) (C 2 ) 




(9) 
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obtain all the matrix elements of vertex operators, and solve the difference 
equation for the simplest case to obtain the two point trace function. 



2 The ABF models 

In this section, after reviewing the corner transfer matrix (CTM) method, we 
introduce the vertex operators (VO's) of the ABF models. We deduce the 
commutation relations among them, express the correlation functions as traces 
of VO's, and derive q-difference equations for the traces. The arguments rely 
on the assumption that the relevant operators are well defined in the infinite 
lattice limit. In the last subsection we relate this heuristic construction to the 
representation theoretical formulation proposed earlier in |l6[ , showing that the 
latter provides a mathematical model for the VO's satisfying all the expected 
properties. 



Definitions 

We consider the ABF models on a square lattice [1], and associate to each site 



j a local state variable Oj G lj = 1, . . . , L — 1, where L is a positive integer 



(> 4). Next, we draw oriented lines on the dual lattice, and associate a spectral 
parameter Qh with each horizontal line and Qv with each vertical line: 



Ch 



Ch 



Ch 



Qv 



Qv 



Qv 



Qv 



To each configuration of local variables (ai, 02, 03, 04) surrounding a face, 
we associate a Boltzmann weight. The latter will be non- vanishing only if the 
configuration is admissible: if the local sites j and j' are connected by an edge 
then <jj — ay = ±1. The Boltzmann weights depend on Q = Qv/Ch, and also 
on a parameter x such that < x < 1, which can be thought of as temperature 
variable: x — > corresponds to the low temperature limit; x — > 1 corresponds 
to the critical limit. We shall restrict our attention to the region 



< x < ( < 1 
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which corresponds to regime III in Q. Since x will remain constant, we will 
display only the dependence on £ and suppress the dependence on x; for instance, 
the CTM below will be written as A^iC). 
Set p = x 2L , and define 



e> p (z) = {z;p) 00 (p/z;p) 00 (p;p)o , 



(z;q 1 ,...,q k ) 00 = (1 - 

ni,...,nfc=0 



With these notations, the Boltzmann weights for the admissible configurations 
are 



0"4 = I4 



&3 = h 



W 



&2 = h 



14 l 3 



C (C = Cv/Ch) 



W 

w 



I l±l 
l±l l±2 



I 

l±l 
I 



l± 1 
Z 

Z±l 



1 c _e p (x 2 c- 2 ) 



K (CV e p (x 2 ) ' 

1 9 p (x 2 'C ±2 ) 

«(C) ©p(^) ' 

1 x 9 p (C 2 )9 p (x 2 (' ±1 )) 

«(0C e p (x 2 )e p (x 2 o • 



(10) 

(ii) 

(12) 



Here the normalizing factor k(C) is so chosen that the partition function per 
site is equal to 1. From the standard inversion trick we find 



«(0 



HO 



e P (x 2 r 2 ) 
e p (x 2 ) ' 



(x 2 C ,P,x 4 )oo(pa; 2 C ,P,% )oo(a; 4 C 2 ;p,a; 4 )oo(K ,P,z )c 



(13) 



(14) 



(x 2 C~ 2 ; p, x 4 )oo ( P x 2 C- 2 ; p, x^Uix^ 2 ; p, x 4 )™ (pC 2 ; p, ^ 4 )c 

Note that «;(£) = ^(xC -1 ). For later use we set 

9i = O p (x 2 '). (15) 



Instead of taking the C,h , Cv 1° be a U the same, we could equally well assign 
spectral parameters that can vary independently from line to line. The resulting 



G 



model is Z-invariant in Baxter's sense. Later we will take advantage of this 
freedom to vary some of the spectral parameters. 

The basic properties of the Boltzmann weights are as follows: 

Initial condition 

l) = huh- (16) 




Unitarity relation 



2>0: / 2 k'M' 4 si<) 



= s, 



h.h- 



(17) 




s h,h 



W 



Crossing symmetry 



^3 h 



c 1 = w 



h h 



9h 



(18) 



Cv 

Yang-Baxter equation 



J2 W 



I h 



h 






h 


H 

-6= 






9h 








9h 


h 




1 


h 










'h 

v h 


h 
I 


Ci/Cs) 



h h 



C1/C2 



'is / 

z 6 /1 



C 2 /C 3 j(i9) 
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\ h I h 

C2 Ci 



Corner transfer matrices 

Let us label the vertices on the lattice as (j'1,.72) (ji,j2 S Z) so that j\ is 
increasing to the left and j'2 is increasing in the upward direction. For i = 0, 1 € 
Z/2Z and m = 1, . . . , L— 2, we denote by C^ m+1 the ground state configuration 
such that 

a Ui ,h) = m+1 if j x + .?2 = i + 1 mod 2, (20) 
= to if ji + j 2 = i mod 2. (21) 

In the following discussion we choose and fix to, and suppress it in our notation. 

Consider the NW-quadrant of the lattice consisting of M + 1 rows and M + 1 
columns, with the vertex (1, 1) at the SE corner: 



1'2 



I'M 



On the boundary shown by bullets, the variables are fixed to have the same 
value as in C^ m+1 . Let us choose a configuration of the local variables 1' = 
. . . , l' M ) (resp. 1 = (Zi, . . . , Im)) on the vertical (resp. horizontal) half line 
(1, 1), . . . , (1, M) (resp. (1, 1), . . ., (M, 1)). Calculate the partition function of 
this quadrant while fixing these variables. The (1, 1') matrix element of the CTM 
A^\Q is defined to be this partition function. Here i refers to the choice of 
the boundary configuration. The matrix element is defined to be zero unless 
li = l[. It is known that in the limit M — > oo the CTM takes the simple form 

where £)W is the CTM Hamiltonian: a £-independent operator whose spectrum 
is contained in the set {0,1,2,...}. The symbol ~ indicates equality up to 
multiplication by a scalar. Such a scalar factor is irrelevant in the computation 
of the correlation functions. 
Set 

= m if k = i + 1 mod 2, 
= m + 1 if k = i mod 2. 

Consider the vector space whose basis elements are admissible configura- 
tions 1 = (li,h,...) satisfying the boundary condition 

l k = f$ if k » 1. 

Let us write |1) when we wish to emphasize that it is a vector in this space. 
Formally the CTM (() or D« act on , and their eigenvectors are certain 
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infinite linear combinations of the |1). In our heuristic approach, we identify TC^ 
with the vector space spanned by these eigenvectors. 

(i) (i) 

We denote by S{ , G\ the operators 

5«|1) = I 1 |1), Gf\\)=g h \\), (22) 

where gi is defined in (|l5|). By construction Z)W commutes with s[ l K Hence 
each eigenspace T-Q 1 ' of with eigenvalue I is invariant under the action of 
£>W, and we have = ®{Hf } . 

The CTMs for the NE, SE and SW quadrants A[ l) ((), A { 2 i} ((), A^(C) are 
defined similarly. Let i?W be the diagonal matrix acting on 

oo 

flWii) = fti?S-|l). (23) 

Then the crossing symmetry implies that 

4 l) (C) ~ RU-(x/0 DW , (24) 

4 l) (C) ~ rW-C^-R®- 1 , (25) 

4 l) (C) ~ G« • (x/0 D(i) ■ R®- 1 . (26) 

Therefore, we have 

A® (C)4 l) (C)4 l) (C)4 4) (C) - Gi° • x 2Dit) ■ (27) 



Vertex operators 

To be able to write down expressions for the correlation functions, we need to 
introduce the VO's 

Notice the difference between the vertex operators introduced above, and those 
defined in the previous section. The point is that the VO's discussed in the 
previous section are defined in the context of vertex models. The VO's that we 
discuss in this section are defined in the context of face models. We will refer to 
the former as VO's of the vertex type and to the latter as VO's of the face type. 
The relation between them will be discussed further below, in the subsection 
" Construction by representation theory" . 

Graphically the (1, 1') matrix element of <j>( l+1 > 1 ) (£) is defined to be the prod- 
uct of Boltzmann weights as follows 
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Ch 



Ch 

h I3 

Ch 

h 1' 2 

Ch 

h\ p[ 

V 

Cv 

The CTM's, their Hamiltonians, and the vertex operators act on specific 

(i) 

subspaces of the full space of states. The action on a certain subspace Hi will 
be indicated by the subscripts I, I', as follows: 



We shall argue that, with an appropriate choice of a scalar factor, the VO's 
enjoy the following basic properties: 

Homogeneity 



^ <i+1) o $(*+m)(c') o c D(z> = & i+1 ^(C/0- 



(28) 



Commutation relations 



h 



C1/C2J $ (i ' i+1) (Cl)/ 4 ,i3* (i+1 ' i) (C2) i3>i2 

<j> (M+1) (c 2 ) Wl <f (4+M) (Ci) W2 . 



(29) 



Normalization 



£< ? * 2 $ (i ' i+1 H<k,* 2 <fr (i+M) (Ck,* 1 



(30) 
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To deduce (28), consider the following: 



B 



m 3 _ 



»3 



m 2 _ 



n 2 



h = P3 ?3 = P2 k = mi = Pi h 



ni 



The figure shows the 

(i+l) 



composition of (£) and The (p, m)-component of A£ (£) is 

represented by the part A and the (m, n)-component of <£>( I+1 ' 4 ) (£) is represented 
by the part B. The composition of A followed by B (with (to 2 , m 3 , . . .) summed 
over) can be viewed as the (1, n)-component of A^(Q, except that the product 

^4 +1) (C)$ (i+M) (C) is acting as H {i) -> H {t+1) while ^(C) as H« -> To 
obtain an equality, and TiS l+1 ^ must be intertwined. In view of the initial 
condition (^) for the Boltzmann weights, the required intertwining action is 
exactly what $( l + 1 '*)(l) does. Therefore, we have 

,(<+1) (Oj a $ (i+1 ' i) (C)i a ,i 1 = • ^ (i) < 



-4 



The normalization (|3 
metry (Q8j) . To motivate 



can be shown using unitarity (|17|) and crossing sym- 
29T) we make use of the Yang-Baxter equation. 





V V 



In the above, 
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the face added in the first figure is pushed up to infinity, using the YB equations, 
and because of the normalization by k its finite effect will disappear in the 
infinite limit. 

Now let us relate the VO's to the correlation functions. As an illustration we 
discuss the simplest situation. Look at the figure below. We wish to calculate 
the local probability of the adjacent variables <j\ and a 2 taking the values l\ 
and l 2 , respectively. We denote this probability by P^ih,^)- 

Divide the whole lattice into 6 pieces. The main parts are the CTM's. The 
remaining pieces give the graphical definition of two types of VO's 

(c) . n {i) H (i+l) $ (M+i) (c) . _^ H {i) 



(t+l,t) 



(0 



A {l+1) 


(0 



















4\o 















<T 2 



(71 

















4\o 



$ (M+D (C) 



The operator $^ +1,l) (C) is the same as we introduced. In 

much the same way as for the CTM $1, M+1) (C) can be written in 

terms of Namely using the crossing symmetry ( |i"8| ) we have 



4> 



(0 



i? W$(- 4 +i)( C ) G ('+i)^+i)-i i 



(31) 



where Gf and i?« are defined in {M> and p 
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Using these operators, we have 



^' <+1) (c)h,^4 i+1) (c) i2 4 +1) (c)/ 2 $^ +1 " ) (c)/ 2 ,h)- (33) 

By ( p8[) and ( |3l| ) we can change the order of the product 

4 M+ %k^4 l+1) (ck4 l+1) (ck 

in d|), to find 

Q {i) {h,h) = 9h9i 2 X 

^(x^U^^xOhm^^HOizm)- (34) 

If we choose the normalization ( |30| ) then the denominator t QW(Zi, Z2) in 
( |32"| ) is equal to tr W (i> ^G^a; 2I?< J J. In a similar way, it is straightforward to 
relate the general correlators to traces of products of VO's. 

The g-difference equation 

Suppose that n is even and the one-dimensional configuration [1%, I2, ■ ■ ■ , l n , l n +i] 
is admissible: lk+i = Ik ± 1, arid l n +\ = l\. We call such a configuration a cyclic 
path of length n and denote it by [h,. . . ,l n ]. Consider the vector space V n 
having the cyclic paths of length n as basis. If li, . . . , l n do not satisfy the above 
admissibility condition we regard . . . , l n ] to be a null vector. On V n we define 
an operator Wk(() (2 < k < n) by 

w k (0[h,...,i n ] = J2 w { l )l 1 i k k +l I <) [h,...M-i,i' k M + i,---M- 



(35) 



Define further 



F, 



,(Cl, - .,Cn) = tr w w (x 2D >U^+ l \a) lul2 ■ ■ ■ ^•HCnhn.l,). 



(36) 



14 



The commutation relation ( |29| ) implies 

'fc+1 

= F h!...dJCl, ■ ■ ■ , (fc+1) Cfc, ■ ■ ■ , Cn)- (37) 
If we set 

FW(Ci,...,Cn) = E f, £.. )lB ^'-'-'^)[ i i."-»U (38) 

il,— ,l n 

then (^) reads as 

W k+1 ((k/(k+l)F^(Cl,---Xn) = F«(Ci,...,Cfc+i, Cfc, ...,Cn). (39) 



Using the cyclic property of the trace and the commutation relations we 
have the difference equation: 



FW(Ci,...,ar j Cfc,...,Cn) = w fc (a_i/a^a)~ 1 ---w a (Ci/^a)~ 1 x 

cw„(c fc /C„) • • • w^ +1 (Cfc/Cfc+i)F (i+1) (Ci, . . . , Cn), (40) 

where C is such that C[/i, . . . , l n ] — [l n , 1%, . . . , £ n _i]. 
Construction by representation theory 

We now proceed to the mathematical construction of VO's. We shall use the 
following notations Q: 

U = Uqi^s^) denotes the quantized affine algebra of type with q = —x 
as the deformation parameter. 

The Chevalley generators are written as e^, /j, U — q hi (i — 0, 1) and q d . 

U' will denote the subalgebra of U generated by a, fi,ti (i = 0, 1). 

V(X) is an irreducible [/-module with highest weight A and highest weight 
vector | A). Ao, Ai denote the fundamental weights. 

V is a two-dimensional representation of U q {sl2) with a standard basis: 
U-}, and V z denotes the affinization of V with the spectral parameter z. 
Note that the grading operator d has a well-defined action on V(A) or V z . Finally 
PR(z) is the i?-matrix, i.e. the intertwiner PR(zi/z<i) ■ V Z1 <£> V Z2 — > V Z2 <S> V Zl , 
which is normalized in such a way that it sends v+ ® v+ to itself. 

In contrast with the situation in the case of the XXZ spin chain, and the 
8-vertex models, here we are interested in VO's of the face-type. The building 
blocks of the latter are still the VO's of the vertex type: 
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$* v (z) : V(X) — > V([i) ® V z . (41) 

Here A,/z are dominant integral weights of level k > 1. Strictly speaking we 
need to complete the spaces properly; but in the sequel we shall not dwell upon 
this point (see 0, Q). A non-trivial VO of vertex type (^|) exists only when 

def 

(tx = X± = A± (Ai — A ), in which case it is unique up to a scalar. We normalize 
them as 

^(z) T \X) = \\ ± )+0(z) 

where the components of ( jil]) are defined by $^ (z) = ^ r $^(z) e g> v e . They 
enjoy the basic properties (fll6|, Appendix 2): 



z d o^(4oz^^(z'/z) E , (42) 



PR(z 1 /z 2 )K v (z 1 )^ x V (z 2 ) 



, \<" v J ipk{pzi/z 2 ) 

J2 ^ (±1+£)/2 $a ± {x- a z)- e ^ (z) £ = 5 Al/ ^ x id y(A) . (44) 



Here we have set p = a; 2 ( fe + 2 ) . 

{x i z;x 4: ,p) QO (z;x 4: ,p) c 



(x 2 z; x A ,p) 



2 



A± = 1 {p r *X 2 ]p) QO 

9x Mp) (p^;p)<» ' 1 j 

and if A = (k — a)Ao + &Ai with Ao, Ai being the fundamental weights, then 

r_ = (a + l)/(k + 2), r+ = 1 — r_. The formulas for w\ are given in |Q, 
cq.(A2.2) wherein q = —a:. 

In the particular case where A has level k = 1, ( pd| ) gives rise to an isomor- 
phism In this case ( |43| ) specializes to 

PR(z 1 /z 2 )^J i (z 1 )^ V (z 2 ) 

— &A-iV( v \if> A *+ lV ( v \ v l-i ( x4z ', x4: )o° ( x2z 1 ', x4 )oo — 
- $ A i+1 ^2j« Ai /_4.-l._4>i r r 2 7 . T 4^ (-2-2i/z 2 ), 

and = 9^ +1 = (a^ 4 ; ^ 4 )oc/(£ 2 ; ^ 4 )oo- Here and till the end of this subsec- 

tion, i is understood to be either or 1. 

Erratum. In the r± in the right hand sides of eqs.(A2.3-4) are to be 
corrected to r-v. 
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Fix L = k+2, £ = (L-m-2)A +(m-l)Ai. Set A; = (L-l-l)A +(l-l)Ai. 
Define 

= Homy (V(Aj), V(0 ® ^(A,)) . (46) 

Namely, H,^ is the space of all U' -linear maps (i.e. those commuting with 
the action of U') f : V(Aj) -> 7(0 <8> V(A,-). Because of the E/'-linearity, / 
is completely determined by the image, say v = f(u\ t ), of the highest weight 
vector u\ [ in V(Aj). Conversely, any vector u in V(£) ® V(Ai), such that it is 
killed by the and has the weight A; (modulo the null root), gives rise to such 
a [/'-linear map. Thus we have the identification 

Jif ={»£ V(0 <8> V(At) | = 0, t,v = q^'^v Vj = 0, 1}. 

The definition of the vertex operators of face type is best described in the 
language of maps, rather than vectors. 
The VO's of face type, 

(i) 

are defined as follows. Given / £ 7i l , consider the composition 
V(Xi) V(£)®V(Ai) 

idm> A . i+lV (z) 

V(0 ® V(A i+1 ) ® 7 



(©^W? +1) (8)V(A,/))®Vi. (47) 



From (^) it follows that there exist unique maps fi> € such that the 

following is true: 

(idy (0 ®*£ +lV (*))o/ = £ (/r®idy)o*£' v (z). (48) 

1<2'<L-1 

We define S^ 1 '^)/-,; by 

* {<+M) (*)jv(/) = /j'. (49) 

It is immediate to see that (z)//^ is nonzero only when Z' = I ± 1 and 

1 < I' < L- 1. 

The following relations are simple consequences of the defining relations ( f48| ) , 
( p9| ) and the corresponding properties (p2|), (E3), (E3) of the VO of vertex type. 

z d Q p+M)^')^ o z -d = §(i+^i)( z '/ z ) VM (50) 
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5 (M+1) (^) i4 ,Z 1 $ (i+1 ' i) (^)! 1) J 3 = E 5(M+1) (^)^^* (i+1 ' i) (^) i3 ,J 5 



W 



zi/z 2 ) , 



(51) 



^ x (lTl)/2 5 A !±1(I ,( M+1)(zx - 2)M±i$(l+M ) (z)i± ^ = ^A ;+1 x . d (g2) 



Here 



W 



h h 



A,: A i+ i 
Ai+i Aj 



1 T77 1 ( ^3 

S(«-V2) w *^A Il A i4 



(53) 



where k — L — 2, p = x 2L , and k(£) is given in (|l4|). From the formulas (A2.2), 
[fi6| one can verify that 



W 



I4 l 3 



= w 



h h 



X Fljl3 F hl2 ( -(l 4 +h-h-h)/2-l 

J F uh Fi lh 



F ll+1 = T p (i), F n _ 1 = r p (i-~) 



(54) 



Here we used ©-(@). 

Now, we return to the ABF models. Let h (h € Cho © Chi © Cd) act on 

h o / — / o h. It is known that the 



H*' by the adjoint action, / S 



principally specialized character of lif*, i.e., tr— (,t)X 2p where p = Aq + Aj, is 

equal to tr ^x 20 ^' . Therefore, we identify H,^ with Wj* and £)W with —p. 
Furthermore, with the identification 



we find that the properties (pgj), (p9h, (p0|) are satisfied by )(C)zz' where 



= {p-x\p)U^^)UP,pZc/{P^x\p)l {x 2 ^ 4 )oo x (1 -p)( L+1 )/ L . 



The Kashiwara-Miwa model 

In this subsection, we wish to give a brief treatment of the Kashiwara-Miwa 
(KM) series as an example of edge- interaction models ^| . The discussion 
of the VO's in the previous subsections can immediately be applied here, since 
an edge-interaction model can be regarded as a special case of a face-interaction 
model. 
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Take an integer N > 2. If TV is even (resp. odd), set N = 2n (resp. 2n+ 1). 
Consider the same square lattice as in the ABF model. For the KM model, we 
consider a local variable <jj taking its values in Z/NZ U {•}. For a E Z/NZ we 
define an integer a* such that < a* < n, a* = a or — a. The admissibility 
condition we impose on (<jj , ay ) is such that one and only one of Oj , ay takes the 
value •. Thus a half of the variables are 'frozen' to the state •. The Boltzmann 
weight is given as follows: 



a 






• 


















• 


\ 


/ 


b 



w 



a • 
• b 



C \a*-b*\ 

,)=^ r MC;a ! 6) 



Cv 



Ch 



W 



• b 
a • 



C =w 



a • 
• b 



x/C 



9b 



where 

9a 

h(£;a,b) 
/(C;«) 

«(C) 
«(0 



1 

Cv 



Q x 2N{-x 4a )/i(l;0,0) 



.9. = 1, 



a .2w(-l)ft(af;O,O) 

/(C;a-6)/(-C;a + 6), 

a*— 1 n — 1 

n e^cr 1 ^ 1 ) n ©,-(c^ +i ) 



„2N . J1N\2 



I— a* 
J2/-2. „4 „4\ 



{x™;x™) 
1 



, X 5 X )oO C 5 *^ 7 *^ )c 



(55) 



Qr 2jV + 2 C 2 ; x 4 , a^ooQr^+y 2 ; x 4 , 



(iV = 2n), 
[N = 2n + l). 



With these Boltzmann weights, we have (|l^-(|l8|) for the KM model. (In 
order to adjust the crossing symmetry to the same formulation as in the ABF 
case, we have modified the parametrization (3.3) in M.) The possible choices 
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of the boundary conditions are similar to (j20|), (pi|): 

= i 1 Ji + h — i + 1 m °d 2, 
= • if ji + j2 = i mod 2. 

(For technical reasons we restrict to the simplest case.) Then all the con- 
sequences of these equations are equally valid in the KM model, including in 
particular the (/-difference equations. Unlike the ABF models, however, we do 
not have a representation theoretic picture of the KM models. In the following 
section, restricting to the Ising case N — 2, we give a mathematical model of 
the VO's in terms of the free fermion algebra. 

3 Vertex operators in the Ising model 

The Ising model is the simplest special case of both the ABF and the KM 
models. In this section we wish to reexamine it in the framework of this paper. 
The special feature about the Ising model is that we are interested in is that 
one can diagonalize the CTM's using Jordan- Wigner fermions, and obtain the 
VO's explicitly. For the reader's convenience we shall repeat the formulation of 
VO's in the context of the Ising model. 

Boltzmann weights 

We choose to work in terms of the edge formulation of the model. In the 
literature, the Boltzmann weights are usually given in terms of the anisotropic 
coupling constants K and L as 





We work in the ferromagnetic low-temperature regime K, L > 0. Following 
[0, let the coupling constants be parametrized as 

sinh(2_ft') = — isn(m), cosh(2i"Q = cn(iu), 
sinh(2L) = ik ns(iu), cosh(2L) = ik~ 1 ds(iu). 

where u — Uv — uh is the difference of (additive) spectral parameters: uy 
attached to the vertical lines and uh attached to the horizontal lines. The 
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sn and cn are the elliptic functions with half periods /, il' corresponding to 
the modulus k — (sinh2i^ sinh2L)~ 1 (||, Chapter 15). We have also used 
Glashier's notation ns(u) = l/sn(u), ds(u) = dn(u)/sn(u), etc. The region we 
consider is then 0<k<l, 0<u<I'. Notice that changing u to I' — u has 
the effect of exchanging K and L (crossing symmetry). 

To compare this with the KM model for N — 2 (see the last figure in Section 



1), let us write w±{Q = W ( a * t),w±(0=W[~ " C)(a,b = 0,1) 



6 



6 



a • 



where ± is chosen according to whether a — b or a ^ b. Then we have 

«(C) 

w-(C) = ^^ x2 ^ 2 ^ 8 )°o(^ 6 C 2 ;x 8 )oo, 
w±(0 = -w±(x/£), 



x. 



X 

(AV 4 ^ 4 )oo(^ 4 r^ 4 ^ 4 )c 

{ x j x )oo 
(x^ 1 , X^)oo 



Identifying x = exp(— ttI'/2I), ( — exp(—nu/2I) we find W-(()/w + (() = 
cn (iu)+isn (iu). Hence, up to the common scalar k(C), the KM parametrization 
agrees with the one for the coupling constants K 7 L. 



Corner transfer matrices 

There are two possible ways to divide the Ising lattice into 4 quadrants. The 
choices depend on whether the site common to all quadrants carries a state 
variable or does not. Thus there are two CTM's to consider, and the complete 
space of states of the model divides into two sectors. This is analogous to the 
situation in the Ising CFT on the annulus, where the space of states divides into 
Neveu-Schwarz (NS) and Ramond (R) sectors, that are half-integrally and inte- 
grally graded, respectively. This will be made explicit below. The possibilities 
are: 

1. The center of the lattice carries a state variable. In this case, the NW- 
quadrant CTM has the form 
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Since the state variable at the center can be either +1 or — 1, this sector fur- 
ther subdivides into two sub-sectors. Once again, this subdivision is analogous 
to the situation in conformal field theory, where the NS sector subdivides into 
two subsectors. 

For that reason, we refer to this sector as the NS sector and write the space of 
states as H NS = Tt NS+ (B'H NS ~ , with the ± referring to the mentioned subsectors. 
The corresponding CTM will be denoted by ^4 S (C) = C 1 ^- Here D NS is the 
CTM Hamiltonian, that the YB equations guarantee to be independent of (. 
The subsectors H NS,± are invariant under the action of D NS . 

2. The center of the lattice is not on the physical lattice. In this case the 
NW-quadrant is 




Once again, by analogy with conformal field theory, we will refer to this 
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sector as H R and the CTM as Af(Q = ( D . 

The CTM's for the rest of the quadrants can be obtained from the above 
by crossing symmetry. To be precise, in the definition of the CTM's we have to 
make a definite choice regarding the boundary conditions at infinity. We choose 
the + boundary conditions. 



Diagonalization of CTM 

The starting point of the diagonalization procedure is to formulate the model 
in terms of Jordan- Wigner (JW) fermions: 

Here the indices of the Pauli matrices refer to the lattice sites labeled by 1, 2, ■ ■ ■ 
from the central site upwards or to the left. They satisfy 

Expressed in terms of the JW fermions, the CTM Hamiltonians are 

j>0 j>l 
j>0 j>l 

Since these are quadratic in the fermions, it is in principle straightforward to 
diagonalize them. Namely we would like to find a set of fermions 

oo oo 
3 = 1 3 = 1 

that diagonalize the adjonit action of the CTM Hamiltonians 

[^ ns ,^s ]= _ 2 ^n Sj [D K^ ] = _2r^. 

As it turns out the eigenvalues are labeled by r £ Z + \ for the NS sector and 
by r £ Z for the R sector. We shall also require that these fermions obey the 
anti-commutation relations 

[<^ S , = Vr S r+sfi , [tf, = VrS r+ s,0, Vr = X 2r + X~ 2r , (56) 

where the normalization factor r\ r is put in for convenience. Conversely the JW 
fermions can be expressed as 



R -6 R . 

r jt'v 
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The determination of a^f and a R is described in ||. Substituting the 
above expansions in the commutation relations and using the expression of the 
Hamiltonians in terms of JW fermions, we obtain second order linear differential 
equations for the generating functions 



(t)=E^'. a R (t) = $> R i'. 

3 = 1 3=1 



For ( p6| ) to make sense we demand that the coefficients a^ s and a R be square 
summable. Changing variables by t — \f~k sn(«) and solving the equations under 
this regularity condition we obtain 

/ IT f 7TT* TTT' \ 

°r S (0 = y 2J ^sc(u) cos(— u) — ivfcsd(u) sin(— v) J , (57) 
a R (t) — sd(v) cos(-^-v) — isc(v) sin(-^-i>)^ . (58) 

They are single- valued and holomorphic for |i| < \/\/~k. 

The Fock spaces 

The 'Hilbert spaces' of states 7i NS and 7i R are the Fock spaces of <p^ s and </> R 
(r ^ 0) over the vacuum states |vac) NS and |vac) R , respectively. The dual Fock 
spaces have dual vacuum states Ns(vac| and R,(vac|, normalized as usual. 

The vacuum states satisfy the usual creation and annihilation conditions 
with respect to cff s and R : 

NS (vac|0, NS = (r<0) ^ s |vac) NS = (r > 0), 

R (vac|^ =0 (r < 0), </> R |vac) R = (r > 0). (59) 

The CTM Hamiltonians are to be normalized so that 

£> NS |vac) NS = D R |vac) R = N s(vac|D NS = R (vac|£> R = 0. 

They give a grading of the spaces 7i , H R by non- negative integers: 

d^'Zi-^.Q.d-.even d€i'Zi>Q,d:odd 



H R = H 

d£'Zi>0,d:even 



R 

d I 



where = {ti£ H NS | L> NS v = dv} and likewise for H R 

Besides the creation/annihilation operators, the fermion algebra in each sec- 
tor contains an extra one ipQ , R : they commute with the CTM's, anticommute 
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with the creation-annihilation fermions, and satisfy (ipQ S ) = 1' (^o ) = 1- 
Hence their action should fix the vacuum vector up to sign. The choice of the 
+ boundary condition corresponds to 



W |vac) NS = |vac) NS , 



Ns(vac|^ s = Ns(vac|, 
R (vac|0o = R (vac|. 



^ R |vac) R = |vac) R , 
In particular H NS,± is the eigenspace of vpQ S with the eigenvalue ±1. 



Vertex operators 



As for CTM's, there are two types of VO's depending on which sector they act 
on: 



$*s V (C):W NS 



n NS . 



Here the superscript V indicates that the VO is 'vertical'. Sometimes the NS- 
subsectors are indicated by a = ±; thus $nL(0 = $ ns (C)| w ns<,, ^r^^OO = 
p*$NSV(Q where pa denotes the projection onto H NSa . 

There are also the 'horizontal' versions of VO's ^nsVCO* $ r S<tH (0- Tnc 
following figures show the (s, s')-elements of these VO's where s = (si, s 2 , ■ ■ ■) 
and so on. The NS-subsectors are indicated by a. 





S\ = G 



R,V 
NS 



NS,V 
R 
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S-2 



Si 



R.H 

NS 




The crossing symmetry implies 

$ NS, CT (0 = $NS,<r(>/C), $R (0 = *r ' 0/0- 

Arguing similarly as in the previous section (the derivation of (pq)), we see that 
the ^-dependence of the VO's is given by 

*S&(0 - c DR ^laK° NS > (60) 

$ NS )CT ,V (C) = r ^ S $NS,^V (1)c D^ (61) 

We will adopt the normalization 

R (vac|<&£s V + (l)|vac) NS = i, NS (vac|$Jf + ' v (l)|vac) R = 1. (62) 



Intertwining Property 

The Boltzmann weights enjoy the initial condition 

ttf±(l) = l, u>+(l) = l, «)_(1) = 0. 

Hence at ( = 1 the VO's become particularly simple: e.g. <£> R S,V (1) carries a 
vertical spin configuration to an identical one on the next left column. However 
it is not the identity because the domain H R and the range 7i NS are different 
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spaces. Using the JW fermions in each sector, the preservation of spins can be 
equivalently stated as the intertwining property 



^NS $ NS^V (1) = $ NS,- CT ,V (1)v ,R +i (63) 

a$f- v (l) = < S ^ V (1)^- (64) 
In the same manner we find the following for the other types of VO's: 



^N + S i$ NS, CT .H (1) 



NS.-ct.H 
R 



«cr$ R (1), 



Along with the normalization condition ([62]), these commutation relations with 
the fermions characterize the VO's uniquely, and all their matrix elements can be 
determined. The working is described in Appendix A. We obtain the following 
formulas for the general matrix elements: 

R (vac|< • • • «s V (C)0 NS fcl ■ • • ^t|v«> OT /C a <- , x-.-.-^+»i+-..+»») 

IK- ii(y~- ' ii v - ii a, n v 



j<j' i,j i<i' 



ii n (-v^T- 2fe *- i/2 7 fcI - 1/2 ) x n 1 1 v fc n 



Here we have set 



(65) 



(66) 



Xk > 1 = 1 _ x 4(k+l) = ~ Xl - k > X -k.-l = X k ,h In = UA. x i\ n n ( 67 ) 
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Vacuum-to-vacuum two point functions 

In what follows, we shall use the shorthand notations 

$ NS )CT) v (c) = $ff(c)) $ R g V (c)=<Ma 

When no confusion may arise, we shall also drop NS, R and write |vac) for 
either |vac) NS or |vac) R , etc. 

Using the formulas for the matrix elements (|65|)-(|66|) one can in principle 
compute the vacuum expectation values of products of VO's. The expansion 
in powers of x suggests that the two point functions are given by the following 
infinite products: setting £ = (1/C2 we have 

^(vac|$ CT (C 1 )^(C 2 )|vac) = { f ' C;a f )o ° F + (C 2 ), (68) 

(vac|$+(C 1 )<i> + (C 2 )|vac) = F_(C 2 ), (69) 

(vac|< / s 2 $-(C 1 )<i>-(C2)^ /2 |vac) = (C 1 + C)F-(C 2 ), (70) 

(vac|<f+(Ci)$-(C 2 )^ /2 |vac) = tx'^F+iC 2 ), (71) 

(vac|< / s 2 $-(C 1 )<i> + (C 2 )|vac) = -ixW^F+i?). (72) 



Here 



F+(z) = 
F_(z) = 



(x Z, X , X )oo 



We shall derive these equations using the formulation of the Ising model as a 
special case of the ABF models (see the end of this section). 
Remark. In the language of the string theory, the vertex operator ^^(z -1 / 2 ), 
for example, is expressed as 



$ns( z 1/2 ) = R( vac l : expF : |vac) 



NS' 



Y =-\ E X m ,nV K m v\z m+n + E X « 



,„R NS m-r 
2 / / --m,»r-mr-n- 1 / y - L m, — r t p~ m L P r z 

m,n>0 m,r>0 



m>0 r>0 r,s>0 

where we have set, for m, n £ Z and r, s £ Z + |, 

For the notation R,(vac| : expY : |vac) NS , see the literature on the fcrmion 
emission vertex, for example 20 1. For x — 1, the vertex operators ^{z^ 1 / 2 ) 
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and $ CT (z- 1 /2) 

are the primary fields with the confomal weight 1/16 of the 
c = 1/2 minimal model, up to some power of z. It can be verified, by noting 
the identity below, that ([38]) and (p^) correspond to the four and two point 
functions of spin fields respectively jcT 

{x x )oo i^x £ ; x ) QO 



1-1/4 

1/2 ]X * J ' 1 + x 2 ^ V 3/2 



1/4 - 1/ '-A^) + ^( 1/ t 3 J A -yye 



where 



n-i 



denotes the basic hypergeometric series. 

Unitarity and commutation relations 

We have the following relations: 

J2f>Ax<W(0=g R xrtnn, (73) 
$ CT «)<i> CT (C) = .g NS xid W Ns,.. (74) 

The reason is that, from the intertwining properties of the VO's, the LHS can be 
shown to commute with the fermions. Since our Fock spaces are irreducible, they 
must then be acting as scalars. The scalars can be determined by considering 
the vacuum-to-vacuum expectation values given above. In this way we get 

R _ l x ; x ; x )oo NS _ l x ? x i x loo 

.gNS tx n n{x™) {-x^U' 

Using a similar reasoning we can derive from (|68|)- ( |72"| ) the following com- 
mutation relations: 

* ff (C 2 )^(Ci) = $ CT (Ci)^(C2)^(Ci/C 2 ), (75) 
QAbWiCi) = X>^(Ci)*^(G0*w(Ci/6)- (76) 

which are what we expect. 
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Difference equation 

Put 

= tr W H (x 2D ($+(Ci)<i> + (C2) ± $-(Ci)$-(C 2 ))) = G*(±G/C 2 ), 
Gf (d/C 2 ) = tr wN s (x 2I3 ($+(C 1 )<i> + (C 2 ) ± $-(Ci)$_(Ca))) - Gf (±d/C 2 ). 

From the commutation relations we then deduce that 

W+ (C)Gf (C) = G^ 2 C), 

G^(x 4 C) = w + (x*Q (w + (() ± «?_(0) G^(C). 

The coefficient in the second equation can be factorized: 



tD+(C)±tD-(C) = 



(ix^£, X^)oc (i^C ^5 *^^)oo 

(x £ , X )oo (x £ , X ; X )oo (»E C , X , X ) 



r«2/ fc -2. ^4^ (^4^-2.-4 ™4\ /'™2/ fa 2. ™4 ™8"\2 
, x ^oo V x S> ) ) *^ /OO V x S> 7*^ 7 X /oc 

Using this we find the following solution: 



(±x 3 c;x 4 ,* 4 )oo ^rvvi. /(o 

" ' " (±xC; x 4 , x- 4 )^ (ix^-i; * 4 , * 4 )oo /(i) ' 



G*(C) = G^ 



/-~2/-2._4 _8 _8-\4 / 10/— 2._4 _8 ^,8^4 /„4/-2._4 „8\ 



/V4/-2-T.4 ~4 -8^2 f~12/*-2. ™,4 ™4 -8^2 ^2/2. ™4 -8^3 ^8^-2.^,2 „81 
^X 1, , X , X , A ^X l, , x , x , x ^ m ^X ,x ,x ; 00 V- 1 ' S j x , x ; c 

The normalization constant C is determined by setting Q = x and using the 
properties ([73|). 



^ /y» 6 ■ .yi4 rt-i4\ 

^X , X , X ^/oo 

Since the solution of a difference equation is determined only up to pseudo- 
constants, it is necessary to verify (e.g. by examining the analyticity) that the 
above solution does give the quantity defined through the trace. We have not 
proved this but checked it to some order in x. 



8 vertex model at the Ising point 

As is well known, the 8 vertex model at a special value of the parameters can 
be regarded as the superposition of two non-interacting Ising models. Let us 
consider a checkerboard lattice whose faces are shaded alternatingly. Place two 
kinds of Ising spins on the faces, one on the shaded faces and the other on the 
unshaded ones. In both cases we assume the -I- boundary condition. Place also 
8 vertex 'spin' variables on the edges taking values ±. Consider a vertical (resp. 
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horizontal) edge separating two faces that carry the spin variables cri, 02. We 
let the edge varibale e have the value 01 02 if the left (resp. upper) face is shaded, 
and — <7i<72 if otherwise. In this way the Ising configurations are mapped to the 
8 vertex configurations with anti-ferroelectric boundary conditon. 

The space of states for this doubled model is the tensor product of the Fock 
spaces H° = H NS <g> H n , Ti 1 = H R <g> H NS . This can be seen by combining two 
Ising configurations to form an 8-vertex configuration: the Ising configurations 
must belong to two different sectors. 

The following operators act on the doubled space H°: ® id, tpQ S ® 
which are mutually anti-commuting. For brevity, we refer to them by <ff,cffi, 
respectively. Similar conventions apply to H 1 . 

For i = 0, 1 e Z/2Z, we define the VO's : W -> W +1 by 

<(C) = £^U0^ S '- ECT ' V (0, 

<J 

*L(0=E*S 8, " v (o®*K^(o. (77) 

(T 

Using the commutation relations for the Ising VO's we find that ( |T7| ) satisfies 

^ + i £2 (C 2 )'i>£ 1 (Ci) = E ^(Ci/C2) £l£2; e;4^ +l£ ;(Ci)'i>5 1 (C2). 

e 1 ,e 2 

The i? matrix is the one given in ]l5[ specialized to g = a; 4 , except that our C -1 
is the C = Csu there; explicitly, in the notation of jTJ| , 

— — = w-(Qw+(0, -77- 4 =w + {Qw-(Q, 

— — = w + (C)w+(C), -77- -r = W-(C)w_(C). 

The matrix elements of VO's can also be obtained in a straightforward way. 
We give below the generating functions of the one- and two-particle matrix 
elements using 

<Hz) = NS (z) + R (z)= E ^z- T . 

Set C= 1, and write $ £ = $S_ e (l) or ^l^ 1 )- 

w, (x 2 z) oa (x 2 w) oa (z~w)(l-x- 1 ^B) 

($+0(z)0M> = Woo Woo (v^WW^Fw < 78 > 

=ix- 1 /»(Vi- VS) ( f ^ ( f T f ) °° , (79) 

(-Sjoo Woo 

Ws (a: 4 ^ -1 )^ (x 2 w) c 
z (x 2 z^ 1 ) oc (w) 

C 



(d>(z)$ + d>(w)) = (1 + ^=) ) °° v ,.„/ , (80) 
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U>(z)S-Mw)) = -ix- l / 2 — {xAz ^ w)o ° (z-w)(x-VHP) 

\[Z {x 2 Z 1 ) OQ (ty)^ [y/z - X\/w){\fz - X 1 \/w) 



<0(«)0(u;)$h 



_ (81) 
1 (x A z~ 1 ) oa (x A w~ 1 ) oc (z - w)(x - yfzw) 

y/zW (x 2 Z~ 1 ) oa (X 2 ^- 1 )^ (y/z ~ Xy/w)(y/z - X~ X %Jw) ' 

(82) 

, n /o \fz — \/w (x^z" 1 )^ (x 4: w^ 1 ) nc 

(ftzMw)*-) = -ix 1 ' ^ 2 _ 1 °° 2 _ 1 °° . 83 

yjzw {x 2 z 1 ) <X3 (aHw 1 ) 0o 



4\ 



Here we used (z)oo = (z;x 
Ising model in the RSOS formulation 

In this subsection we examine the Ising model formulated as an ABF model 
with L — 4. Using the representation theoretical construction we shall calculate 
the two point functions for the VOs. 

For our purpose it is useful to study first the symmetry under the Dynkin 
diagram automorphism. Let v denote the algebra automorphism of U' given 
by v(ei) = ei_i, v(fi) = fi-u v{U) = ii-i- For A = (k - l)A + IA X set 
A = IAq + [k— l)A\ . Then there exists a unique isomorphism (which we denote 
also by v) of vector spaces V(X) — > V(A) such that f(|A)) = |A) and v{x\u)) = 
v{x)v(\u)) for all x € U' and \u) £ V(X). To make manifest the symmetry under 
v let us introduce the VO of vertex type in the principal picture: 

$^ y (C) = c (A1_A ' p) x (id®c p )5A y (C 2 ) 

where V (z) is the VO @. In our case p = A\ — Aq. The power of £ is supplied 
in order to retain the normalization < I , ^ ±1 V (C)|A;) = |A;±i) <8> w T + 0(C)- We 
have 

Given / G Home// (V(A), V(£) <8> Vfa)) we write v(f) = (^yjo/o^ 1 G 
Hom^ (V(A), V(£) (8> F(f?)) . Define the VO of face type in the principal picture 

^(C) by 

(id v( a ® ^ V (C)) o / = E ® id) ° $r(0- (84) 

A' 

Then we have the symmetry 

^ ) (C)(/) = ^(^' n) (C)K 1 (/)))- 

Now we focus attention to the Ising model. Taking k — 2,we write "* (C) 
for 3> A A A'f iA< '^(0 (recall that A; = (2 — £)A + £Ai). Exhibiting the j dependence 
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explicitly we set 

n[ r ' l) = {v £ V(Aj) <g> V(Ai) I e s v = 0, t„v = q {h °' Xl) v (s = 0, 1)}, 
so that 

V(A 3 ) ® V(Ai) = ® V{Xt). 

l=i-\-j mod 2 

Clearly z/ induces the isomorphism 7Yp'^ c± H^_/ }1 ^ ■ We shall identify 

h (0;0) =w n S + j W (0;0) =H N S _ 5 ^1) ^(1:0) 

and set 

|vac) NS - |Ao) ® |A ) e ^ 0;0) , 
|vac) Rli = |A 4 )®|Ai_ i )eW| iil ~ <) . 

Defining the dual vectors in a similar manner, we wish to compute the two point 
functions of (£)• From their defining property (84) it follows that 

<Ao|$^(Ci)$t y (C 2 )|Ao) 

= Ns(vac|$S ;01) (Ci)$lo ;10) (C2)|vac) NS x (X \^f (Ci)*£ v (C 2 )l A > (85) 
(A^i^CCi)^, 1 -^^)^,) = £ a,i- i (vac|Sg- i;il - i HCi)^ i ( 1 1 " i;1 "" ) (C2)|vac) R)1 _ i 

/=0,2 

x(A 1 |$^^(Ci)$^(C 2 )|A 1 ). (86) 

On the other hand, the following can be obtained by solving the qKZ equa- 
tion for VOs of vertex type. 



(x 4 CV 4 )c 



(T&CI.^A rf&\ f^lifl.^A rf%\ 



(£ 8 C 2 ;^,z 4 ) c 



<Axi*r (<o*r (gom = w ,;j ) , / 14 " ;c4 ^ >< 

1/4 V - 1/4 ; x 8 , ^/f i * 10 C 2 ) *, 

Here £ = Ca/Cii an d the sign ± is chosen according as i — 0, 1 or I = 0, 2, 
respectively. Using the symmetry 

R^vaclSE' 01 ^)^ 
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we can solve (|§), (||) for the two point functions of the $$ H>) (0- With thc 
identificaiton 

*THo = * + (c 1 ), $r i) (c) = $-(r 1 ), 



we then recover the formulas given in (|6 
Creation and annihilation operators 

To conclude this section we shall briefly discuss the creation and annihilation 
operators for the row-to-row Hamiltonian. We shall introduce two kinds of such 
operators, and discuss the relation between them. 

So far we have chosen to work with the + boundary condition for the CTM. 
To take the excited states into account, it is necessary to consider both of the 
two boundary conditions cry = ± (j 3> 1). We denote the space for the CTM 
with those boundary conditions by Tt± respectively, where Ti. refers to either Ti R 
or 7i NS . The two spaces are identified by i — Y\( >0 erf (in the naive picture), 
which is the disorder operator at the origin. The action of the VO's on 7Y_ is 
defined by $°"(C)|w_ = i° 3>~ CT (C) ° <- and $ ff (C)|w_ = to $_ ff (C) o l. 

A priori there are 2x2 possible choices of the boundary conditions for the full 
space on which the row-to-row Hamiltonian acts. For simplicity we shall fix the 
boundary condition for the left half to be +, and set J r +,± = Hom(H±, 7i+) ~ 
H+ <g> H*±. The dual space is ± = Hom(H+,7i:±) ~ H% O H± where the 
coupling between J-+ : ± and T*L ± is chosen to be (f\g) = tr^ ± (fg) for / S J 7 ! ± 
and g € 3-+,±- In this picture the ground state vector and its dual are |vac) = 
(vac | = x D . 

The first way to introduce the creation and annihilation operators is via the 
explicit diagonalization of the Hamiltonian using fermions. Define 



We find that the creation and annihilation operators ip(z), (p*{z) : J~±,+ —* ^~^,+ 
are given by 

tp(z)(v) = ip(xz) o v o l — v o ip(zjx) o l, 

ip* (z)(v) = ip(z/x) o v o l + v o ip{xz) O L, (87) 

where v G J~±.+- They satisfy 

(p(z)(\vskc)) = 0, [<p(zi),<p*(z2)]+ = S(z 1 /z 2 ). 

The second way is to utilize the transfer matrices Tj} s (C), t^ s (Q in the 
representation theoretical picture, as is done for the XXZ Hamiltonian and 
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other models H,Ir|, Let us set 

t*5(0(«) = (O^E $ (cv o«o$(cr, 

(7 

^ S (o(«) = cO _1 E *(cr ° « ° *(0*. 

If we define 

<fi(z)(v) = <j){z/x) O V O L, (f)*(z)(v) — <j){xz) O V O L 

then we have the following commutaion relations: 

t* s (04> NS (z) = f(zC 2 /x 2 )$ R (z)r§ s (0, 

Here f(z) signifies the function defined in Appendix A, (p8|). 

The two types operators denned above are related, for example, as follows. 

0*(z)(|vac>) = <£*(z)(|vac>) 

0*(*i)(0*(* 2 )(|toc») = ^(z 1 )(^(z 2 )(|vac))) - ,5(x 2 z 1 /z 2 )|vac). 



Acknowiedgement. The authors would like to thank Fedor Smirnov for inspir- 
ing discussions. One of the authors (M.J.) is grateful to Department of Mathe- 
matics, ANU and Brian Davies for kind invitation and hospitality. This work is 
partially supported by the Grant-in- Aid for Scientific Research on Priority Ar- 
eas, the Ministry of Education, Science and Culture, Japan, and the Australian 
Research Council (ARC). 

A Matrix elements of the Ising VO 



In this appendix, we shall outline the derivation of the matrix elements of VO's. 

R 



Let us consider $^ S ' V (1). Define f(z) and c r by 



/(*) = -VkMv) = -^- 1/2 t^¥t£^ ■ 

where v is related to z by z = exp(niv / 1) . Using the Fourier expansion f(z) 
J2 r e-+Z °r z ~ T on \ z \ — 1 an d the formulas (|57|), (|58| ) it can be shown that 

oo 

E] a ^/ a -s j + l^Vs 1 = c r-s- 
3=1 
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With the aid of this we rewrite the intertwining property (163) as follows 



NS $ NS,V (1) = £ Cr ^NS,V (1)0 K 

In terms of the generating functions 

NS (z) = J2 tf**^, (89) 
reZ+i 

R (z) = $>v. ( 9 °) 

seZ 

we have NS (z)$ R S ' V (l) = /(z)$r S ' V (1)<?!) R (z), or equivalently 

NS (z)^ v (l) = ci>f' V (l)0 R (z). (91) 

Here we have introduced auxiliary operators 

NS (z) = /_(z)0 NS (z)= £ ^ s z-™, 

«eZ 

R (z) = / + (z)0 R (z) = X! (9 2 ) 

neZ 

according to the factorization f(z) = f + (z)/f—(z), with 

= , f_( Z ) = ix -l/2 z l/2 & iz ~ 1 > 3Si )°° 



(x^Z,X^)oo (x^ Z ^,3^)oo 

Notice that we have 

NS (vac|^ s = [n< 0), 

R |vac) R = (n > 0), ^|vac) R = |vac) R . (93) 

In view of ( [E)l| ) no confusion may arise if we abbreviate expressions like 

NS <vac|0 NS (z)d>f ' V (l)|vac) R = NS <vac|< S ' V (l)^ R (z)|vac) R 

to Ns( v ac|0(z)|vac) R and so on. From the above properties, it follows that 

(z — w) (z + w — (1 + x~ 2 )zw) 



NS 



(vac|0(z)0(ui)|vac) 



(z — x 2 w)(z — x 2 w) 



Note that this formula contains Ns(vac|</>(z)|vac) R = 1 as a special case w = 0. 
The general matrix elements are given by Wick's theorem. 
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The other case $^g V (l) can be treated quite similarly. The intertwining 
property becomes 

where f(z) is the same as (|88j). Defininig auxiliary operators 

R '(z) = /i(z)0 R (z), ? NS '(z) = f' + (z)<f> NS (z) 

with ^ 2 14 

/' (z) = ix 1/2 z~ 1/2 ^ Z,X /' (z) = — 2 

we have 

^'(z)^ V (l) = $« s V (l)0 NS '(z), 

R (vac|</> R ' = (n < 0), R (vac|</> R ' = R (vac|, 0^ s '|vac) NS = (n > 0). 

We hnd 

(z — w)(z + w — 1 — ir 2 ) 



(vac|0'(z)^'(w)|vac) NS 



(z — a; 2 w)(z — x 2 w) 



The matrix elements relative to the creation/annihilation operators in each 
sector can be derived from the above generating functions. Let X mn and j m be 



as in (67). Using the formula 



( z . x 4) = 1^ ?»* 
[Z,X Joo n=Q 

we obtain 

Ns(vac|$^ s ' v (l)0 R m (/) R I Jvac) R = X n ^ n j m j n (m,n > 0,m + n > 0) 
NS (vac|^ s +1/2 $f ' V (l)0 R n |vac) R = X_ m _ lA ^ 2m+1 / 2 7 m 7„ (m,n > 0) 
NS (vac|0^ s +1/2 ^ i/2(f) NS,v (1)|vac)R = Xrn+1/2 , n+1/2 x 2m+2n+1 lmln . (m,n> 0) 

Likewise we have 

R (vac|0 R R $ R s V (l)|vac) NS = X m , n x 2m+2n lm ln (m, n > 0, m + n > 0) 
R (vac|0 R $ R s V (l)0 N s_ 1/2 |vac) NS = X m ,_ n _ ll2 ix 2m - X l 2 ^ n (m,n> 0) 

R< v ac|$NS V ( 1 ) ( / , -L-l/2^-n-l/2l vaC >NS = X m+l/2,n+\/2 1 )7m7n • (m, 71 > 0) 

Finally, by making use of the identity (N even) 



Pfaffian 



n f 



1 XiXj J 1 < i j< N x<i<j<JV ^ 

we arrive at the general expressions (pa), (pq). 
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B Spin correlation functions 



The spin correlation functions for the Ising model are given as traces of VO's. 
They can be calculated using the intertwining property and the unitarity of 
the VO's alone, without referring to the formulas for the matrix elements. In 
this appendix, we illustrate this by the simplest case of the nearest diagonal 
order-order and disorder-disorder correlations. 

Let us introduce some notations. Setting z — exp(iriv / 1) wc define 



/o NS (*) = v^-(«) = vS £ ^ v < 

/ R (z )= y^dn(,)=y^5>-v. 

V 7T V 1 *-L 

They are related to f(z) = —y/ksn(v) given in ( |8S| ) by 

/ N V^)/(z) = if*(z), tf(x- 2 z)f(z) = -z/ NS (z). 
Define further 



In terms of them and the generating functions (|89|), (|90[), the intertwining prop- 
erties of the VO's are summarized as follows: 

</> R (z)<MC) = f(zC 2 )<s>-AC)^ s (z), 

*$°(Q = ^(C)V'f(C), 
Rewriting the properties for the horizontal VO's we have further 

$ ff K)0 R (z) = /(zC 2 )0 NS W<f- CT «), 
$ ff «)^ NS (z) = /(zC 2 )0 R M<i>-.«), 

<rW(x£) = z^ S (C)<f^(C), 
a<Z> a (x() = ^(0<M0- 

The order-order and disorder-disorder correlators for the nearest diagonal 
neighbors are given respectively by 



(94) 
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£ CT tr W R (x w * « 2 W (Ci (C2 ) 

3 NS ff R tr H H(x 2CR ) 



(95) 



Here the ^ are the spectral parameters attached to the intermediate lines as 
shown below. 

Consider (|94|). Using the intertwining properties and the unitarity of the 
VO's we have 

^^ ff Ki)$ CT (Ci) = E $ -«i) $ff (^)V' R (Ci)i = 5 R V' R (Ci)i- 

cr cr 

Hence, setting ( = C1/C2, we find 

^aa'$ ff '« 2 )^Ki)<i> ff (Ci)^(C2) 

(J(j' 

= ff R ^a'$ ff 'K 2 )^ R (Ci)i^(C 2 ) 

= 5 R //f-/o R (^)E^ <7 '(^)^ R wci 2 )^(c 2 ) 

J Z7TIZ * — ' 
cr' 

= ff R /^-/o R w/(^/c 2 )^ s (^/c 1 2 )E^^'^)^(c 2 ) 

cr' 

- -W* 3 / ^/ o R (,)/(./c 2 )0 NS (./c 2 )^ NS (c 2 ). 

To take the trace we invoke the following simple lemma: 

^ — ^ L=S*(x 2 z 1 /z 2 ), 6*(z) = J2z r , 

tr « E ( r ) r tz 

tr„ N s fa; 2DNS f > NS (z 1 ) f > NS (z 2 ) N ) 

Here the <5 NS (z) contains half odd integer powers in z, but in the course of the 
computation only integer powers appear. With the aid of these formulas, ( p4[ ) 
becomes 

The case of can be treated similarly. We find that it is given by 

2r 

'Ir <» 
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